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Abstract 

Considering the quintom model with arbitrary potential, it is shown that 
there always exists a solution which evolves from u) > — 1 region to uj < — 1 
region. The problem is restricted to the slowly varying potentials, i.e. the 
slow-roll approximation. It is seen that the rate of this phase transition 
, only depends on the energy density of matter at transition time, which 

■ itself is equal to the kinetic part of quintom energy density at that time. 

\Q ' The perturbative solutions of the fields are also obtained. 

o 

oo ■ 1 Introduction 

O 

o 



Recent astrophysical data indicate that the expansion of the universe is accel- 
erating [1]. This acceleration can be related to the presence of a perfect fluid 
with negative pressure, known as dark energy, which constitutes two third of 
our present universe. 



£~ I ■ One candidate for dark energy is the cosmological constant: a constant quan- 

tum vacuum energy density which fills the space homogeneously, corresponding 
to a fluid with a constant equation of state, EOS, parameter u>d = —1 (in 
this paper u>d and ui denote the EOS parameters of dark energy component 
■ and the universe, respectively). This model suffers from fine tuning and coin- 

cidence problem [2]. Alternatively, dynamical homogeneous fields, e.g. scalar 
fields with canonical and non-canonical kinetic terms and with various poten- 
tials, have been proposed as the origin of dark energy [3,4]. In these models, 
in contrast to the cosmological constant model, the EOS parameter may vary 
with time. The EOS parameter satisfying — 1 < uj < —1/3, which determines 
the quintessence era, can be achieved by introducing a normal scalar field </>, 
known as quintessence scalar field [3]. 

Due to the fact that the observations have shown some mild preference for 
an equation of state parameter u < — 1, the phantom scalar field a with a wrong 
sign kinetic term was introduced in the literature [4]. Depending on phantom 
scalar field potential, different solutions such as asymptotic de Sitter, big rip, 
etc. may be obtained [5]. 
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Some astrophysical data seem to slightly favor an evolving dark energy EOS 
parameter u>d and show a recent Ud = — 1 crossing [8]. As the dark energy model 
with a quintessence field has an equation of state u>d > — 1 and for the phantom 
ghost field we always have uod < — 1, the phase transition from quintessence to 
phantom era does not occur in models with only one scalar field, either phantom 
or quintessence [6]. To study such a transition, we need models which are at 
least composed of two scalar fields [9], known as hybrid models. One of these 
models is the quintom model which assumes that the cosmological fluid, besides 
the matter and radiation, is composed of a quintessence and a phantom scalar 
fields [7]. In [10], a phase-space analysis of a spatially flat Friedman-Robertson- 
Walker universe containing a barotropic fluid and phantom-scalar fields with 
exponential potentials has been presented. This model has late time phantom 
attractor solution. In [11] the same calculation has been done by introducing 
an interaction term between the phantom and the quintessence fields. The 
quantum stability of quintom models is also an important problem whose some 
aspects have been discussed in [12]. 

Recently a new view of quintom model, known as hessence model, has been 
introduced. In this model, the dark energy is described by a single field with 
an internal degree of freedom rather than two independent real scalar fields. 
Hessence model can avoid the difficulty of the Q-ball formation which gives 
trouble to the spintessence model [13]. The evolution of u in this model has 
been studied in [14] and [15], via phase-space analysis. Although this model 
allows u> to cross —1, but it avoids the late time singularity or the "big rip". 

In the recent paper [16], we have investigated the necessary conditions for 
occurrence of the u = — 1 crossing in the quintom model ( which also implies 
the Ud = — 1 crossing ), in special the transition from quintessence to phantom 
phase. This transition has been checked for two specific quintom potentials, the 
power-law and exponential potentials, and has been shown that the u = — 1 
crossing really exists in these examples. It has been assumed that the fields 
have slowly varying behavior, i.e. the problem has been solved in slow-roll (SR) 
approximation. 

In this paper, we are going to prove that this phase transition occurs for an 
arbitrary quintom potential, as long as we consider the SR approximation. We 
think that this is an important result which is in agreement with the present 
data. 

The scheme of the paper is as follows. In section 2, we briefly review the 
main results of [16] and write down the necessary conditions needed for occur- 
rence of the quintessence to phantom phase transition in the quintom model. 
In section 3, by solving the Friedman equations for arbitrary potential around 
the transition point and in the SR approximation, it is shown that the desired 
transition occurs for the general quintom model. The general form of SR con- 
ditions, needed for consistency of the equations, is stated in terms of the time 
derivatives of quintom fields. In section 4, we obtain the explicit perturbative 
solution of quintom fields for arbitrary potential, from which the SR approxima- 
tion can be expressed in terms of the derivatives of the potential, which is more 
physical. This is discussed in appendix A. We check our perturbative results 
with the exact expressions known for special cases. 

We use the units h = c = G = 1 throughout the paper. 
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2 Transition conditions in quintom model 



Consider a spatially flat Friedman-Robertson- Walker universe with scale factor 
a(t), filled with (dark) matter and quintom dark energy. The evolution equation 
of matter density p m is 

p m + 3Hj m p m = 0, (1) 

in which j m = 1 + u> m . w m is the equation of state parameter of matter field 
defined by uo m — P m /p m , where P m is the pressure of matter field. H(t) = 
a(t)/a(t) is the Hubble parameter and "dot" represents the time derivative. The 
quintom dark energy consists of a normal scalar field 0, i.e. the quintessence 
field, and a negative kinetic energy scalar field a, called the phantom field. The 
energy density po and pressure Pd of the homogenous quintom dark energy 
are [7, 10] 

Vfaa), (2) 
and their evolution equations are 

2 i oiri i dVfaa) 

4> + 3H 4> -\ — — = 0, 

dip 

* + 3Ha- dV ^ =0. (3) 

OCT 

The Friedman equations, obtained from Einstein equations, are 
tt2 871" . 
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Pd = 






2 ( 





4tt • 

— {4> 2 -& 2 +2V(^c7) + 2p m ], (4) 



and 



H = -4tt( Pd + p m +P D + P m ) 

= -Att(4> 2 -<7 2 + J m p m ). (5) 

Note that eqs.©, Q and, (O are not independent. The equation of state 
parameter lo = (Pd + P m )/(pD + pm) can be expressed in terms of Hubble 
parameter as 

—!-§£• W 

For H < 0, the system is in the quintessence phase lo > —1, and when H > 0, 
it is in the phantom phase with to < — 1. So crossing the u = — 1 line is, in 
principle, possible in quintom model. 

If we are interested in situation in which the quintessence to phantom phase 
occurs in some instant of time t = to, then H(t) must have a local minimum at 
that time, i.e. H(to) =0. So at t < to, H < and to > — 1 and at t > to, H > 
with lo < — 1. If we restrict ourselves to t — <o << ftj , where ho = H(to) and 
hy 1 is of order of the age of the universe, H(t) can be taken as 

H(t) =h + hi(t - t Q ) a + 0((t- t Q ) a+1 ) . (7) 
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a > 2 is the order of the first non-vanishing derivative of H(t) at t — to and 
hi = ±H( a \t ). H^(t ) is the n-th derivative of H(t) at t = t . The desired 
phase transition occurs when a is an even positive integer and hi > 0. 

For arbitrary quintom potential V((f>, a), we want to study if this situation 
( even a and positive hi ) exists or not. As usual, we restrict ourselves to SR 
approximation in which the first terms of equations Q are negligible: 

4> « H(j), 

a « Ha. (8) 



3 Seeking for transition solution 

Following [16], we are going to find any solution to cqs.(J3])-((5|), when H(t) given 
by eq. (|7j) and the fields vary slowly, i.e. eq.©, with the desired property, that 
is h\ > and a = even. 

Let us begin with eq.® and expand it near to = 0. One finds 

ahi^- 1 + ■■■ = -4tt[/3(0) + $(0)t + •••], (9) 

in which 

(3{t) = j> 2 ~ a 2 + lmPm . (10) 

In terms of the dimensionless variables $ = 0//io, £ = v/ho, R m = p m /hQ, 
t = hot, and Hi — hi/h^ , eq.® can be written as 

d Ft 

aHiT*- 1 + ■■■ = -An[B(0) + — (0)r + •■•], (11) 

dr 

for r << 1 ( t « /iq 1 ). B(t) is defined through 

B(t) = & -t 2 + lm R m . (12) 

We will continue our study with eq. Q , reminding that shifting to dimensionless 
variables is always possible. 

The matter density p m (t) can be found by solving eq.([T]) with H (t) given by 
eq.©. The result is 

Pm(t) = p m (0)e- 3 ^' iot+ ^ tQ+1 >. (13) 
For a > 2, the first equation obtained from eq.© is 

0(0) = 2 (O) - a 2 (0) + 7 mPm(0) = 0. (14) 
For second relation, we first note that 

0(0) = 2(00 - aa) - 3h 0l 2 mPm (0). (15) 
But in SR approximation ([8|), we have 

- (H<j)- 4>)<j)+ (Ha - a)a ~ H{a 2 - <j> 2 ). (16) 

So 

- &a) « h (& 2 - 2 ) o = h Q j m p m (0), (17) 
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where in the last equality we use eq. (|14[) . Therefore, excluding special j m << 1 
cases, one finds 

-4tt/3(0) = 127r^o7™Pm(0), (18) 

which is a non-zero positive quantity. This shows that for all quintom potentials, 
one has 

a = 2 

h X = 67T/lo7mPm(0) > 0, (19) 

which proves the uj = — 1 crossing for all quintom models in SR approximation, 
as long as the remaining relations which are obtained from eq.((3]), are satisfied, 
up to the lowest order, consistently. 

The important observation is that in all quintom models in SR approxima- 
tion, the rate of the phase transition (hi) only depends on the (dark) matter 
energy density p m (0), which itself is determined by the kinetic part of quintom 
energy density feq. (fl"4")l ). 

It is worth noting that if there is no dark matter field, the eq. (fT8|) , and 
therefore eq.JTH]), can not be used directly. In this case, we can not neglect the 
first two-terms of eq. (|TB|) and therefore hi becomes 

hi = -4:ir(<j)4> - <7«t)o. (20) 

The above equation shows that the lu = — 1 crossing is in principle possible in 
dark-matter free quintom models. 

To prove the consistency of the remaining relations, we expand eq.(U]) near 
t = 0, which results in 

hi + 2/ l0 M a + • • • = [(5(0) + 5(0)t + l'5(0)t 2 + •••], (21) 

in which 

6(t)=<j> 2 -a 2 + 2V(0,<j)+2 Pm . (22) 
The first two relations obtained from eq. ([2"Tj) are 

= h 2 , (23) 

which determines ho in terms of V(0) and p m (0), and 5(0) — 0, which can be 
written as 

+ V) t =o = 0. (24) 

3H(a 2 -^ 2 ), (25) 

in which the equation of motion Q in SR approximation has been used. So 
eq. ([24|) can be written as 

U{4> - 3H<j>) - &{a - 3H&) - 3h olmPm ) - 0, (26) 

V / 4=0 

which is reduced to eq. (|14p in SR approximation. 



y(0) + (l- JjL)p TO ( ) 



era - 3/i 7 



But 



V - —<])+—& 
d(h da 
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The last equation that must be checked in the lowest order approximation, 
is one obtained from t 2 -term in the right-hand-side of eq. pTj) . Using eas. (fT5|) 
and and H(0) = 0, it becomes 

47T ■■ • •• .... 

— [<f> 2 - 6h <jxf) - a 2 + Gh aa + $</>—&'&' + 9/i 7 m p m ] t=0 

47T .... 

~ —[4,(4, - 6h Q 0) +&{-&■ + 6h Z) + 9h 2 ol 2 mPm ] t=0 , (27) 

where the SR approximation has been used in the second equality. But the 
natural extension of SR approximation ((8j) to higher order terms is 



d n (f> 




d"- 1 ^ 




<< 




~dF 


d"cr 




d n ~ l a 




<< 


k °dt"-l 


~dF 



(28) 

So the eq. ([27|) reduces to 

4irh Q [2{aa - 44) + 3/io7mPm]t=o = -47r^ 0/ 9(0) = 2/i /u, (29) 

where eas.(|15|l and ©, with a = 2, have been used. Now this is exactly equal 
to the second term in the left-hand-side of eq. pTj) for a = 2. 

Let us briefly express what we have done. We have shown that the eqs.©- 
([5|) have a transition solution when H(t) behaves as ([7]) with a — 2. It has been 
shown that eqs.(j4]) and §5$) result three independent relations: Eq. (fT4|) . which 
relates p m (0) to the kinetic energy of quintom fields, and eqs. flT))) and ([2"5|) 
which determine the Hubble parameters hg and hi in terms of matter field and 
potential. For t « h^ 1 and in SR approximation, it has been shown that all re- 
maining relations are consistent with the mentioned three ones. This completes 
the proof of u) = — 1 crossing of all quintom models in SR approximation. □ 

It is worth noting that the validity of eq. (|28|1 can be verified by following 
argument. Consider the equation of motion © near t = 0. Differentiating these 
equations with respect to t, and using H(Q) = 0, one finds 

d dV 

^ + 3 ^ + d7^ = °< 

d dV 

<T + 3H(j — — — — — = 0. (30) 



dt d, 



a 



Now ignoring the <fi ( a ) term in eq.Q in comparison with Hip ( Ha ), is 
equivalent to ignoring the cf> ( a ) term in eq. (|30|) in comparison with H<p ( Ha 
). Repeating this differentiation procedure justifies eq. ([28|) as the general SR 
conditions. 



4 Solution of equations of motion 

In special cases, like ones considered in [16], one can, in principle, solve the 
equations of motion ^ in SR approximation and studies their physical behav- 
iors. For arbitrary potential, as expected, we can not solve these equations 
without knowing the functional form of the potential, but instead we try to find 
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a perturbative solution, from which one can deduce some physics behind the 
problem. Here we focus on the equation of motion of the quintessence field. 
The solution of phantom field can be simply obtained from final relations by 
replacing the derivatives from <f> to a and transforming V — > — V. 

Consider the equation of motion of quintessence field in SR approximation 
4> « H(j) with H = h + h x t 2 : 



3(/i + M 2 )</> = - 



dV(cf>,a) 



Using the Lerchphi function a, 6) defined as 

®(z,a,b) = J2 

it can be easily shown 



— ' (b + n) a ' 

n=0 v ' 



cl 

df 



Eq.(|3l]l then leads to 



2/in 



felt 2 1 
' ho ' '2 



M 2 ' 



3G(</>) + A, 



where 



G(4>) = - 



(31) 



(32) 



(33) 



(34) 



(35) 



dV/d^ 

and A is the constant of integration. Note that in this stage, we assume that 
dV/dcj) does not depend on a. Expanding both sides of eq. (J33J) around t = up 
to order i 3 , results in 



hi 
3/i 



2 t 3 + 0{t 5 ) = 3[G(0) + G{Q)t + ^G(0)t 2 + ^G{0)t 3 + O{t 4 )}+ A. Ci(ii 



1 

3? 



which leads to 



3G(0) + A = 0, 
1 

h ' 

G(0) = 0, 



3G(0) = 



(37) 



^G(O) 



A. 

3^' 



Using 



where 



A -\jal' B ={w)«- C =W)„ (39) 

one can obtain 0(0), 0(0), and 0(0) from eg. (l37j) . from which 0(i) is obtained 
up to order t 3 as following 

«*> " * ( »» - + ^ - 4 - 2Al ) '* + 0<t,) - m 

This is our desired expression. 

As an example, we consider the exponential potential 

V = wie Al * + v 2 e X2 ° . (41) 

Eq. (|40|) then results in: 



, wiAie Al ^° u 2 A 3 e 2Al *> 2 wiAie Al *° /2u 2 A?e 2Al *> „, \ , A 

^ - ^-^k^ + ^k§-' — w — 2/11 J * + ° ( * 

(42) 

But it can be easily shown that the above expression is same as one obtained 
from expanding the exact expression: [16] 

*(*) = 77 T > ^ ' ^ ' 



Ai ^ ^ A? f**(=fe^, 1, |) + 2dfc, 
in which 

ci = -^e- x ^°. (44) 

At the end, let us discuss a possible ambiguity. It may be argued that as we 
take H(t) up to order t 2 in eq. (f3"Tj) . keeping the £ 3 -terms for 0(t) is not correct. 
To answer this question, we keep the i 3 -term in H(t): 

H(t) = h + h x t 2 + h 2 t 3 + ■■■ (45) 

In SR approximation, we have 

^- = 3G(0) + A. (46) 
Using (|45p . the left-hand-side of eq.JH]), up to order t 3 , is 

A_I^3_ 1^4 + ... m 

h 3h 2 Ah 2 + 1 l) 

So the coefficient h 2 only contributes to i 4 -term and therefore the expansion 
(|40| is correct. 
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Appendix A: general SR condition for potential 



In this appendix, we want to obtain the analogous of SR conditions (|28[) for 
quintom potential, which is probably more physical. This can be done by using 
the explicit solution (1401) . For simplicity, we assume that the parameters A, B, 
and C, defined in eq . (f3"9")) , are positive real numbers. Using (|4"0"|) . we first note 
that the condition \<p\ << ho\<p\ simply results in 

B=(?£) «K. (A.l) 



This constraint can be also obtained by noting that in SR approximation, 
we have 3ft <$(0) + d/di(W/<90)| o = ( see eq.^Sfl} ), which results 0(0) = 
-{l/3h )(d 2 V/d(j) 2 )4>\o. Then \4>\ « h \4>\ reduces to ea.lTOl. 

The next conditions can be obtained by nothing that the consistency of t 2 - 
term of eq. (|2"Tj) can be achieved if the {<j> 4>)o ( and (& '&')$ ) term in eq. (|2"Tj) is 
negligible in comparison with hlp m (0) term. But eq. (fT4"l) shows that p m (0) ~ 
2 (O), so the equations are consistent if 

hj> 0) O «0 2 (O). (A.2) 
"o 

Using eq. (|40|) . eq. (|A.2p reduces to 

1 -A -A [B 2 +AC , \ A 2 
2hi « 



h 2 3h 3h 2 \ 9h J 9k 2 , 



1 f B\ 2 hi 1 A C , /A oN 



Now eq. p^)) implies 



hi pm{0) 4> 2 (0) 

7,3 ~ U2 ~ u2 ~ 1.4 ■ 

"o "o "o "o 



(A.4) 



where eas. p^)) and (|4"01) have been used. So eq. (|A.3|) is satisfied if, besides 
condition (|A. 1|) . one has A « /iq and C << /iq. The generalization is obvious: 
the general SR conditions, in the language of quintom potential, are: 

'd n V" 



« K, 



ii 



d n V 
da 11 



« K, (A.5) 



o 

which indicate the slowly varying potential. 
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